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ABSTRACT 



Most models of cosmic chemical evolution predict that the mass-weighted 
mean interstellar metallicity of galaxies should rise with time from a low value 
~ 1/30 solar at 2; ~ 3 to a nearly solar value at z = 0. In the absence of any 
selection effects, the damped Lyman-alpha absorbers (DLAs) in quasar spectra 
are expected to show such a rise in global metallicity. However, it has been 
difficult to determine whether or not DLAs show this effect, primarily because 
of the very small number of DLA metallicity measurements at low redshifts. 
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In an attempt to put tighter constraints on the low-redshift end of the DLA 
metalUcity-rcdshift relation, we have observed Zn II and Cr II lines in four DLAs 
at 0.09 < z < 0.52, using the Space Telescope Imaging Spectrograph (STIS) 
onboard the Hubble Space Telescope (HST). These observations have provided 
the first constraints on Zn abundances in DLAs with z < 0.4. In all the three 
DLAs for which our observations offer meaningful constraints on the metallicity, 
the data suggest that the metaUicities are much lower than the solar value. These 
results are consistent with recent imaging studies indicating that these DLAs may 
be associated with dwarf or low surface brightness galaxies. 

We combine our results with higher redshift data from the literature to esti- 
mate the global mean metallicity-redshift relation for DLAs. We find that the 
global mean metallicity shows at most a slow increase with decreasing redshift. 
For the redshift range 0.09 < z < 3.90, the slope of the exponential fit to the 
binned A^(II I)-weighted mean Zn metallicity vs. redshift relation is —0.18 ±0.06 
counting Zn limits as detections, —0.22 ± 0.08 counting Zn limits as zeros, and 
—0.23 ± 0.06 using constraints on metallicity from other elements in cases of Zn 
limits. The corresponding estimates of the z — Q intercept of the metaUicity- 
redshift relation are -0.74 ± 0.15, -0.75 ± 0.18, and -0.71 ± 0.13, respectively. 
Roughly similar results are obained if survival analysis or an unbinned A'^(HI)- 
weighted nonlinear approach is used. Thus, the A^(H I)-weighted mean metal- 
licity of DLAs does not appear to rise up to solar or near-solar values at ^ = 0. 
This weak evolution could be explained by the fact that our absorption-selected 
sample seems to be dominated by dwarf or low surface brightness galaxies. This 
suggests that current DLA samples, especially those at low redshifts, could be 
biased against more enriched galaxies because the latter may cause higher dust 
obscuration of the background quasars. 

Subject headings: quasars: absorption lines; galaxies: evolution; galaxies: abun- 
dances; cosmology: observations 

1. Introduction 

A great deal of progress has been made in recent years in the observations of high- 
redshift galaxies. The average star formation history of the Universe has been estimated 
from emission properties of galaxies detected in deep imaging and redshift surveys (e.g., 
Lilly et al. 1996; Madau et al. 1996). This emission history of galaxies is tied intimately 
to the history of gas consumption and metal production in galaxies, because the global 
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mean densities of gas, metals, and stars are coupled through conservation-type relations. 
Constraining the history of metal production thus puts constraints on the histories of star 
formation and gas consumption (e.g., Pei & Fall 1995). 

Most models of cosmic chemical evolution, ranging from analytical calculations to nu- 
merical simulations, predict the mass- weighted mean global metallicity in galaxies, given by 
^metals /^gas^) ^ise from a low value at high redshift to a near-solar value at zero redshift 
(e.g., Pei & Fall 1995; Malaney & Chaboyer 1996; Pci, Fall, & Hauser 1999; Somerville, 
Primack, & Fabcr 2001; Tissera et al. 2001). [Here f^^^ais ^^'^ ^gas^ refer, respectively, 
to the mean comoving density of metals and the total mean comoving density of gas in 
the interstellar matter of galaxies, both measured in units of the present critical density 
(pc = ^Hq/StiG).] The present-day mass- weighted mean metallicity of local galaxies is also 
observed to be close to solar (e.g., Kulkarni & Fall 2002). It is therefore of great interest to 
determine whether element abundance measurements in galaxies actually show the predicted 
rise in mean metallicity with time. 

Abundance measurements in quasar absorption systems can directly probe the history 
of metal production in galaxies. Particularly useful in this respect are the damped Lyman- 
a (DLA) absorbers (log A^hi > 20.3). Under reasonable assumptions, DLAs may contain 
enough neutral gas to form a large fraction of the stars visible today (Wolfe ct al. 1995). 
DLAs are our principal source of information about the chemical content of interstellar 
matter in high-2; galaxies (e.g., Meyer & York 1992; Lu ct al. 1996; Kulkarni ct al. 1996, 
1997, 1999; Pcttini ct al. 1994, 1997, 1999; Prochaska & Wolfe 1999; Prochaska et al. 2001a, 
2003a, 2003b; and references therein). 

For a variety of reasons, the abundance of Zn is often taken as the metallicity indicator 
in DLAs. Zn tracks Fc closely, to within ~ ±0.1 dcx, in Galactic halo and disk stars for 
metalhcities [Fc/H] > -3 (e.g., Mishcnina ct al. 2002; Cayrel et al. 2004; Nisscn et al. 
2004). Zn is also relatively undepleted on interstellar dust grains, so that gas-phase Zn 
abundances can give metallicity estimates relatively free of depletion effects. Finally, Zn 
II has two absorption hues at AA2026, 2062 that, in DLAs, often lie outside the Lyman-a 
forest and are usually unsaturated. There are, however, some uncertainties related to the 
nucleosynthesis of Zn and the potential role of ionization corrections, and we return to these 
issues in section 5.1. 

The absorption hues of DLAs sample the interstellar matter of distant galaxies. This 
sampling is completely independent of many properties of the galaxies, including their lumi- 
nosities, sizes, and morphologies. It can be shown (e.g., Lanzetta, Wolfe, & Turnshek 1995; 
Kulkarni & Fall 2002) that, in the absence of selection effects, the quantity f^m^ais/^gl^ 
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galaxies is equal to the A^ni-weighted mean metallicity Z in a sample of DLAs, given by 

This quantity Z may therefore be expected to show a rise from a low value at high redshift 
to a near-solar value at low redshift. 

There has been great debate about this issue in the past few years, with most studies 
advocating no evolution in the global mean metallicity Z inferred from DLAs (e.g., Pettini et 
al. 1997, 1999; Prochaska & Wolfe 1999, 2000; Vladilo et al. 2000; Savaglio 2001; Prochaska, 
Gawiser, & Wolfe 2001b). Kulkarni & Fall (2002) examined this issue further by applying 
various statistical analyses to the Zn data from the literature. Based on a compilation of 57 
Zn measurements (36 detections and 21 limits) from the hterature at redshifts OA < z < 3.4, 
Kulkarni & Fall (2002) found the slope of the metaUicity-redshift relation to be —0.26 ±0.10, 
consistent at f=:i 2 — 3(7 level with both the predicted rates of evolution (-0.25 to -0.61), and 
with no evolution. Recently, Prochaska et al. (2003b) have reached similar conclusions 
(slope —0.25 ± 0.07) for 0.5 < z < 4.7, using Zn measurements in 11 DLAs and Fc, Si, S, 
or O measurements in 110 other DLAs. As in the case of Zn, the elements Fe, Si, S, and 
O also have both advantages (e.g., better theoretical understanding of nucleosynthesis, case 
of detection) and disadvantages [e.g., dust depletion (for Fc, Si), lack of two or more weak 
lines outside the Lyman-alpha forest (for Si, S, O), potential ionization effects (for Si, S), 
potential enhancement with respect to the Fe-group elements due to alpha enrichment (for 
Si, S, O), etc.]. In any case, the results of Prochaska et al. (2003b), based on these other 
elements, are similar to those of Kulkarni & Fall (2002), as discussed further in section 4. 

The issue of metallicity evolution in DLAs is still not fully understood. The main reason 
for the debate about this issue is the small number of measurements available, especially at 
z < 1.5. The samples of Kulkarni & Fall (2002) and Prochaska et al. (2003b) both include 
only 1 measurement each in the redshift range z < 0.6, which corresponds to > 40% of the age 
of the Universe. The low-redshift systems require space-based ultraviolet (UV) measurements 
to access the H I lines at ^ < 1.6 and the Zn 11 lines at z < 0.6. With the goal of putting better 
constraints on the low-redshift end of the DLA metallicity-redshift relation, we have obtained 
UV spectra of three quasars with four intervening DLAs at 0.09 < z < 0.52 using the Space 
Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Here we 
report our new observations of the Zn II and Cr II absorption in these low-redshift DLAs and 
their imphcations for the metaUicity-redshift relation for DLAs. We have chosen to focus on 
Zn as the primary metallicity indicator, rather than combine measurements of Zn with those 
of Si or S, primarily because (a) we prefer to use the more uniform Zn sample, relatively 
free of assumptions about dust depletion, and (b) STIS is much more sensitive at the longer 
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wavelengths of the Zn II hnes than at the wavelengths of the S II and Si II lines for our DLAs. 

Section 2 describes our observations and data reduction. Section 3 describes the results 
for the individual DLAs. Section 4 presents the implications of our results, combined with 
higher redshift data from the literature, for the metallicity-redshift relation of DLAs. Section 
5 discusses the advantages and disadvantages of using Zn, the morphological information for 
the galaxies giving rise to the DLAs in our sample, and potential selection effects. Finally, 
section 6 provides a summary of our results. 



2. Observations and Data Reduction 
2.1. Observations 

Our sample consists of 4 confirmed DLAs with 0.09 < z < 0.52 and logiVni > 20.3 in 
the sightlines toward the quasars Q0738+313, Q0827+243, and Q0952+179. Table 1 lists 
the key properties of these DLAs known from previous studies. The hnes of interest are Zn 
II A A 2026.136, 2062.664, and Cr II AA 2056.254, 2062.234, 2066.161, which occur in the 
range of 2211-3145 A for the DLAs in our sample. 

Our observations were obtained during February- November 2003 using HSTSTIS. Table 
2 summarizes the details of our observations. The STIS near-ultraviolet Multi- Anode Mi- 
croChannel Array (nuv-MAMA) was used for Q0738+313 and Q0952+179. For Q0827+243, 
the CCD was used because of its higher efficiency at the longer wavelengths required for this 
object. The first-order gratings G230M and G230MB were used to get the moderately high 
spectral resolution and wavelength coverage necessary for the abundance determinations. 
The 52" X 0.2" slit was used for the MAMA observations to achieve a good compromise 
between spectral resolution and throughput. For the CCD, the 52x0. 2E1 pseudo-aperture 
(with the default target placement near row 893, ~ 5" from the top of the CCD) was used 
to mitigate the effects of the charge-transfer inefficiency. 

All targets were acquired with the F28x50LP aperture using the STIS CCD in the 
acquisition mode. Acquisition exposures of 3-5 s were used to achieve a signal-to-noise ratio 
(SNR) of > 40, taking care not to saturate the detector. The science exposures were dithered 
using the "STIS- ALONG-SLIT" pattern in 3, 4, or 5 steps of 0.4 " each, so as to help with 

flatficlding and removal of hot pixels. For the CCD observations, the exposures at each 
dithered setting were divided into 3 equal subexposures with CR-SPL1T=3 to facilitate the 
rejection of cosmic rays. An automatic wavelength calibration exposure at the respective 
grating setting was included in each orbit. 
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2.2. Data Reduction 

All data reduction tasks were performed independently by two members of our team, 
using multiple techniques as discussed below. The spectra were reduced using the Image 
Reduction and Analysis Facility (IRAF) version 2.12.1 and the CALSTIS package in the 
Space Telescope Science Data Analysis System (STSDAS) software version 3.1. The auto- 
matic lamp exposures obtained during our observations were used for wavelength calibration. 
For all other cahbrations, the most recent reference files available in the HST data archives 
were used. These include image reference files for darks, pixel-to-pixel flat flelds, low-order 
flat flelds, tables for aperture description, aperture throughput, bad pixels, dispersion co- 
efficients, incidence angle correction, and image distortion correction for both the MAMA 
and the CCD. In addition, the CCD reductions required the bias calibration files from the 
HST archives. For the MAMA, other calibration files included tables for linearity, offset cor- 
rection, scattered light correction, near-UV dark correction, and time-dependent sensitivity 
correction. 

After the initial reduction of each 2-dimensional image, the spectral extraction was 
carried out using the STSDAS task XID. In order to obtain optimum signal-to-noise ratios in 
the 1-dimensional spectra, a variety of extraction heights were tested. Extraction heights of 3 
or 4 pixels were found to encompass all of the recoverable flux without introducing additional 
noise from the background. The results of XID were checked by repeating the spectral 
extractions with the IRAF task APALL. In these latter reductions, the sky background 
level was fltted in the regions on either side of the quasar spectrum, and this flt used to 
subtract the background level. In most cases, second-order Chebyshev polynomials were 
found to provide adequate flts to the traces of the spectra. The results of these various 
reductions were found to be in very close agreement. Finally, for each target, the individual 
1-dimensional spectra at each wavelength setting were combined together. The combined 
spectra were then continuum-fltted using Chebyshev or cubic spline functions of orders 2 to 
6, depending on the shape of the continuum and the presence or absence of broad emission 
lines in the wavelength range observed. Figs. 1-5 show the continuum-normalized intensity 
versus vacuum wavelength for all the targets. The spectral dispersions are 0.09 A pixel^^ for 
the NUV-MAMA spectra, and 0.15 A pixel^^ for the CCD spectra. The resulting spectral 
resolutions are in the range of ~ 21 — 24 km s~^ for the NUV-MAMA data and ~ 29 km s~^ 
for the CCD data. 

Table 3 lists the features measured in the flnal spectra. The central wavelengths and 
equivalent widths were measured using the IRAF task SPLOT. The 1 a uncertainties in 
equivalent widths due to photon noise were estimated assuming Poisson statistics as = 
AA-^iVpja,/ SNR where AA is the spectral dispersion, Np^^ is the number of pixels occupied 
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by the feature, and SNR is the signal-to-noise ratio per pixel in the continuum near the 
feature. 

Another uncertainty often not considered in spectroscopic studies of QSO absorption 
line systems is the uncertainty in the equivalent width (and hence column density) associated 
with the continuum placement (e.g., Savage & Sembach 1991). While the continua of QSOs 
are in general much simpler than the stellar continua encountered in studies of the Galactic 
interstellar matter, there is still some uncertainty associated with the exact placement of 
the continuum. In order to assess this uncertainty, we slid the continuum up and down 
by 0.3 times the 1 a noise level and measured the equivalent widths for these extrema 
(see Sembach & Savage 1992 for a justification of the scale factor 0.3.) . Table 4 lists the 
continuum placement uncertainties in the equivalent widths cr^ thus calculated for the DLA 
metal features. Also listed are the total uncertainties cr^ calculated by adding the photon 
noise uncertainties and continuum uncertainties in quadrature. Further details on the line 
measurements are provided in section 3. 



2.3. Estimation of Column Densities 

Column densities for Zn II and Cr II were estimated by fitting the observed absorp- 
tion line profiles using the program FITS6P (Welty, Hobbs, & York 1991) that evolved 
from the code used by Vidal-Madjar et al. (1977). FITS6P minimizes the between 
the data and the theoretical Voigt profiles convolved with the instrumental profile. The 
atomic data were adopted from Bergeson & Lawler (1993) for Zn II A A2026, 2062 and Cr II 
AA2056,2062,2066, and from Bergeson, Mullman, & Lawler (1994) for Fe II AA2249,2260. 
The oscillator strengths used for these lines were 0.489, 0.256 for Zn II AA2026,2062; 0.105, 
0.078, 0.0515 for Cr II A A2056, 2062, 2066; and 0.00182, 0.00244 for Fe II A 2249, 2260. 

In an attempt to constrain the column densities and effective Doppler b parameters, 
we performed simultaneous single-component fits to the detected lines of Zn II, Cr II, and 
Fe II (where available), varying the column densities N, the Doppler parameter b, and the 
central velocity v together, but assuming a common b and v value for Zn II, Cr II, and the 
weak Fe II fines. We also performed simultaneous single-component fits to both detected 
and undetected lines. The uncertainties in N, b, and v were determined simultaneously in 
these fits. 

In the cases where a 2026 A feature was potentially present, we investigated its contri- 
bution by carrying out profile fits with the Mg I A 2026 line included. To do this, the Mg I 
column density was estimated from profile fitting of the Mg I A 2853 feature in data obtained 
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at the Multiple Mirror Telescope (MMT), if available (Khare et al. 2004). The MMT data 
were obtained as part of a larger survey of clement abundances in DLAs at z < 1.5. These 
data have lower spectral resolution (~ 75 — 95 km s~^) but higher SNR (~ 50 — 90). The Mg 
I A2026 line was found to make insignificant contribution to the observed 2026 A feature. 
Indeed, there is no evidence of absorption at the expected positions of Mg I A 2026 in any 
of the 3 DLAs where the relevant wavelength region was covered. 

The best-fit effective b values were found to be in the range of 11-15 km s^^. However, 
due to the moderate resolution and SNR of our data and the weakness of the lines, it was 
difficult to derive accurate estimates for the effective b values. Therefore, in order to explore 
the dependence of the estimated column densities on the b values, we also performed single- 
component fits to the individual detected lines for several (fixed) choices of b, i.e., 5, 10, 15, 
and 20 km s~^. In general, relatively little difference in column density was found for the fits 
with 6- values of 5-20 km s"^, indicating that the lines would be near the linear part of the 
curve of growth for that range in effective b. To search for any broad, shallow components 
that could be hidden in the noise, we also examined all the spectra after smoothing them by 
a factor of 3. No broad and shallow components could be discerned. 

The only other available information relevant to the effective b values for the DLAs in 
our sample is from the H I 21 cm absorption lines (Chcngalur & Kanekar 1999; Kanekar 
& Chengalur 2001) and/or from lines of Ca 11 or Mg 11 seen in our MMT data for these 
systems (Khare ct al. 2004). While the latter data cannot be used to estimate the 6- values 
of individual components accurately, they provide upper limits to the effective 6- values, and 
indicate a spread in velocities of several tens of km s^^. Further details on the b values are 
provided in discussions of the individual systems in section 3. The total column densities 
obtained with the present resolution are expected to be close to the actual values if very 
strongly saturated components are not present (Jenkins 1986). It seems unlikely that any 
of our Zn II, Cr II, or Fe II lines consist of single components with 6-values less than 5 km 
s"-*^. Indeed, b values < 5 km s~^ for the dominant component contributing to the Zn II 
lines are found in none of the eight DLAs at z < 1.5 and only f» 20% of the 24 DLAs at 
z < 3.4 for which profile fits from high-resolution data are available in the literature (Lu 
et al. 1995; Meyer et al. 1995; Prochaska & Wolfe 1997; Lopez et al. 1999; Pettini et al. 
1999, 2000; de la Varga et al. 2000; Molaro et al. 2001; Ellison & Lopez 2001; Ledoux et al. 
2002; Levshakov et al. 2002; Peroux et al. 2002; Ledoux et al. 2003; Dessauges-Zavadsky 
et al. 2003, 2004). Thus, it would be unusual, although not impossible, in a sample of four 
low-redshift DLAs, to find two or more with b < 5. [Of course, it is possible that the low 
incidence of systems with b < 5 km s~^ is a reflection of the limited resolution (?a 6-8 km s~^ 
at best) achieved so far for DLA spectra with current instrumentation. In principle, if DLAs 
often contain significant narrow, unresolved components, then many existing column density 
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estimates for heavy elements in DLAs could be too low.] It would be desirable to have higher 
resolution observations in the future to understand the detailed component structure. This 
would need a large amount of observing time on a space-based telescope. 

In order to evaluate the effect of continuum placement uncertainties on our column 
density estimates, we slid the continuum up and down by 0.3 times the 1 a noise level and 
measured the column densities for these cxtrema by refitting the lines with FITS6P. This 
was then used to calculate the uncertainty in the column densities due to the continuum 
placement uncertainty, which was then added in quadrature with the errors associated with 
the photon noise, uncertainty in 6-values, and uncertainty in component velocities to derive 
the final errors. 

Finally, we also estimated the column densities using the apparent optical depth (AOD) 
method (see, e.g.. Savage & Scmbach 1991). The column densities obtained using this 
method were consistent with those obtained from profile fitting within the mutual 1 a un- 
certainties. However, the AOD results are expected to be less accurate than the profile-fitting 
approach for our data because of the relatively low SNR and resolution. This is because the 
AOD method weights all points across the observed profile equally, i.e., it does not take into 
account the fact that the observed profile "should" resemble the instrumental profile if the 
lines are unresolved or marginally resolved. Thus, for example, outlying downward noise 
points could spuriously increase the computed AOD column density. We thus adopted the 
profile-fitting estimates of the Zn II, Cr II, and Fc II column densities and combined these 
with the known H I column densities to estimate the clement abundances. For reference, we 
adopted the solar abundances recommended by Lodders (2003). 



3. Results 
3.1. Q0738-I-313 

This quasar {zem — 0.635) has two strong DLAs dX z — 0.0912 and z — 0.2212 with 
log ATjjj = 21.18l[]:[J^ and log A^hi = 20.90;[|;|J^, respectively (Rao & Turnshek 1998). The 
z — ^fl2Vl system was known to be an Mg II absorber and identified as a DLA on the 
basis of EST UV spectra (Rao h Turnshek 1998). The z — 0.0912 absorber was discovered 

serendipitously in the same EST UV data. The quasar also shows 21-cm absorption at the 
redshifts of both the DLAs (Lane et al. 1998; Chengalur & Kanekar 1999). The Zn II and 
Cr II lines of both systems lie outside the Ly-a forest. The wavelength coverage achieved 
with three separate grating settings included all five lines of Zn II and Cr II for each DLA. 
The "unidentified" features near 2155, 2187, and 2217 A are hkely to be C IV A 1548 hues 
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ai z — 0.392, 0.413, and 0.432, respectively, with the corresponding A1551 features lost in 
the noise. 



3.1.1. The z = 0.0912 DLA 

Fig. 6 shows the velocity profiles of the Fe II, Zn II, and Cr II data for this DLA. The 
zero point of the velocity scale refers to the absorption redshift z — 0.0912 inferred from the 
21-cm and optical data. Counting both photon noise and continuum uncertainties, our data 
show weak ^ 2.8 o", 2.7 a, and 3.3 a features, respectively, at the wavelengths of the Zn II 
A2026, Cr II A2056, and Fe II A2260 lines. No other lines of these species were detected in 
our HST spectra. The 21-cm absorption in this object has been observed to be dominated 
by a single component at z = 0.09118±0. 00001 with a full width at half maximum (FWHM) 
of ~ 5 km s^^, indicating a Doppler parameter 621cm = 3.0 km s~^ for cold, neutral hydrogen 
(Chengalur & Kanekar 1999). ^. Based on our complementary MMT spectra (Khare et al. 
2004), we derive effective 6 < 42 km s~^ for Ca II. 

Wc first carried out profile fits using only the potentially detected lines Zn II A 2026, Cr 
II A 2056, and Fe II A2260. The fits were then repeated by fitting these lines together with 
the undetected lines Zn II A 2062; Cr II AA2062, 2066; and Fe II A2249. To be conservative, 
we adopt the results from these latter fits, and hst the derived Zn II, Cr II, and Fe II column 
densities, and the corresponding abundances in Table 5. The error bars on column densities 
and abundances for these fits include uncertainties from Poisson statistics and uncertainties 
in b value, central velocity, and continuum placement. The formal estimate of the best fit Zn 
II column density in this DLA has a relatively large total error, therefore we list it only as a 3 
a upper limit. To explore the dependence of the estimated column densities on the b values, 
we also performed single-component fits to the individual detected lines for several (fixed) 
choices of 6, i.e., 5, 10, 15, and 20 km s~^, which are also listed in Table 5. The column 
densities for these single line fits with b fixed are higher than that for the multiple-line fit 
with N, b, V varied, because the former included fits to only detected lines, while the latter 
included fits to both detected and undetected lines. 

The dashed curves in Fig. 6 show the profiles for the best-fit N, b, and v parameters 
obtained by varying these parameters simultaneosuly for the multiple lines stated above. 
The estimate of b value obtained in this way is 11.2 ± 6.7 km s~^, where the uncertainties in 



^The Zn II, Cr II, and Fe II lines can arise in warmer H I gas or even H II regions, and could thus have 
effective 6-values different from &2icm- (If H II regions are present, the derived abundances [X/H] would, 
strictly speaking, be upper limits.) 
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fe were calculated simultaneously with the uncertainties in N and v. 

For pictorial illustration, we created simulations of the Zn II A 2026 line for solar, 32 
% solar, and 10 % solar metaUicities ([Zn/H] = 0.0, -0.5, and -1.0), assuming the best-fit 
b and v values, b — 11.2 km s~^ and v — 9.9 km s~^. These simulations are shown in Fig. 
7, overlaid on top of our data, and support the profile fitting results that the metallicity 
appears to be low. 

3.1.2. The z = 0.2212 DL A 

This DLA has also been detected in 21-cm absorption. Kanekar, Ghosh, & Chengalur 
(2001) have reported three components in the 21-cm absorption line, the dominant one 
having FWHM = 18.6 km s~^ and comprising ~ 87% of the H I column density. This 
indicates a Doppler parameter &2icm = 11-2 km for most of the cold, neutral hydrogen. 
Our MMT data suggest effective 6 < 38 km s~^ for Mg I and Mg II in this DLA (Khare et 
al. 2004). 

Fig. 8 shows the velocity profiles of the Zn II and Cr II data for this DLA. The zero 
point of the velocity scale refers to the absorption redshift z = 0.2212 inferred from the 
21-cm and optical data. Our data show pa 2.8 a and 3.5 a features at the wavelengths of the 
Cr II A 2062 and Zn II A 2062 lines. These lines lie in the C IV broad emission line of the 
quasar, and therefore have higher continuum SNR than that for the Zn II A2026 and Cr II 
A2056 lines which are not detected. 

We first carried out profile fits using the potentially detected lines Zn II A 2062 and Cr 
II A 2062 together with the undetected line Cr II A2066. The undetected Zn II A 2026 and 
Cr II A2056 lines were not included in the fits because they lie in regions of the spectrum 
with lower continuum signal level. The estimate of the b value obtained in this way is 
12.2 ± 7.7 km s~^. In Table 6, we list the derived Zn II and Cr II column densities, and the 
corresponding abundances. The dashed curves in Fig. 8 show the profiles for these best-fit 
N, b, and v parameters. Also listed in Table 6 are results of single-component fits to the 
individual detected lines for several (fixed) choices of b, i.e., 5, 10, 15, and 20 km s~^. The 
column densities of Cr II for these single line fits with b fixed are higher than that for the 
multiple-line fit with N, b, v varied, because the former included only the detected lines (Cr 
II A2062, Zn II A2062), while the latter also included the undetected Cr II A2066 line. 

We next checked whether the Zn II and Cr II column densities thus estimated are 
consistent with the lack of significant features at the positions of the Zn II A2026 and Cr 
II A2056 lines in the noisier parts of our spectra. To do this, we generated simulated Zn II 
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A2026 and Cr II A2056 lines with the best-fit parameters estimated from the Cr II A A 2062, 
2066 and Zn II A 2062 hnes. At our resolution, the Zn II A2026 and Cr II A2056 features 
would be 7 pixels wide, but are lost in the noise. The 3 a limits on the rest-frame equivalent 
width over these 7-pixel wide simulated features are about 89 mA and 75 mA, respectively, 
at the positions of the Zn II A 2026 and Cr II A 2056 lines. For the best-fit h value of 12.2 
km s~^, this implies 3 a upper limits of log A^znii < 12.81, i.e. [Zn/H] < —0.72, and log 
Ncxii < 13.37, i.e. [Cr/H] < —1.18. These limits are consistent with the values obtained 
from the data for Zn II A2062 and Cr II A A2062, 2066. Thus the Zn abundance in this 
system is likely to be close to -0.70, the value obained from the A2062 line. Nevertheless, to 
be conservative, we adopt an upper limit [Zn/H] < —0.70 for this system. 

For pictorial illustration, we show in Fig. 9 simulations of the Zn II A 2062 fine for solar, 
32% solar, and 10% solar metaUicities ([Zn/H] = 0.0, -0.5, and -1.0), assuming the best-fit 
h and v values, h — 12.2 km s~^ and ^; = 15.3 km s"^. The parameters of the Cr II A 2062 
line are held fixed in these simulations, as determined from the best- fit obtained by fitting 
together the Cr II A A 2062, 2066, Zn II A 2062 hnes. These simulations support the profile 
fitting results that the metallicity appears to be low. 

3.2. Q0827+243 

This quasar {zem = 0.939) has a DLA at ^ = 0.5247 with log A^hi = 20.301^;^^ (Rao 
& Turnshek 2000). The system was known from ground-based studies as a strong Mg II 
absorber at z = 0.5247, and found to be a DLA in HST UV spectra (Rao & Turnshek 2000). 
The quasar also shows a relatively broad and possibly multi-component 21-cm absorption 
feature at ^ = 0.52476 ± 0.00005, with FWHM ~ 50 km s'^ , indicating 621cm ^ 30 km s"^ 
in the cold neutral hydrogen (Kanekar & Chengalur 2001). Our MMT data show Fe II and 
Mn II absorption at z = 0.5249, with an effective 6 < 35 km s~^ (Khare et al. 2004). 

Unfortunately, this quasar was in a significantly fainter state at the time of our HST 
observations than expected from previous UV data. Therefore the SNR in our HST data 
was poorer than expected. The wavelength region covered all five lines of Zn II and Cr II, 
but no hnes could be detected. Fig. 10 shows the velocity profiles for this DLA. The zero 
point of the velocity scale refers to the absorption redshift z — 0.5247 inferred from the 
21-cm and optical data. Although weak features are seen at -20 km s~^ in the Cr II A 2062 
line and possibly at -40 km s~^ in the Zn II A2062 line, they are comparable to the noise. 
Furthermore, there is no evidence for these features in the stronger Zn II A 2026 or Cr II 
A 2056 lines. The 3 a rest frame equivalent width limits in our HST data at the expected 
positions of the Zn II A 2026 and Cr II A 2056 lines are 96 mA and 90 mA, respectively. 
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The corresponding 3 a upper limits on the Zn II and Cr II column densities and abundances 
arc listed in Table 7 for 6 = 5, 10, 15, and 20 km s~^. The differences between the limits for 
6 = 15 and 20 are small, so the limits for higher 6-values are expected to be close to those for 
6 = 20 km s^^. For example, for h = 30, the b value of the 21-cm absorption line, the limits 
are about 0.1 dex lower than for b = 10. The dashed curves in Fig. 10 show the profiles 
corresponding to the 3 a upper limits for b = 10. 

The low SNR of our STIS data for Q0827+243 prevents us from putting tight constraints 
on the Zn and Cr abundances in the DLA absorber along this sightline. We note that based 
on MMT spectra of this target, the abundance of Fe in this DLA is [Fe/H] = —1.02 ± 0.05 
(Khare et al. 2004). The total (gas + solid phase) metallicity could of course be higher. But 
the DLA toward this quasar has a much lower A^hi than that of the other three DLAs in 
our sample. Therefore, unless the metallicity of this DLA is substantially supersolar (which 
would imply a very large amount of dust depletion), it would not significantly affect our 
estimates of the global metallicity of DLAs at low redshifts, which anyway treat this DLA 
as an upper limit. 



3.3. Q0952+179 

This quasar {zem = 1-472) has a strong DLA at z = 0.2378 with log A^hi = 21.321°;°^ 
(Rao & Turnshek 2000). The absorber was known to be a strong Mg II system and found 
to be a DLA in HST UV spectra (Rao & Turnshek 2000). The quasar also shows 21-cm 
absorption at the same redshift (Kanekar & Chengalur 2001). 

Our wavelength setting covered the Zn II A 2062 and Cr II A A 2056, 2062, 2066 lines. 
A separate setting to cover only the Zn II A 2026 line was not obtained in view of the large 
number of orbits needed for this faint quasar to cover even a single grating setting to a 
reasonable SNR. Tables 3 and 4 list the lines measured in the final spectrum. Fig. 11 shows 
the velocity profiles for the Zn 11 and Cr 11 lines in this absorber. The zero point of the 
velocity scale refers to the absorption redshift z = 0.2378 inferred from the 21-cm and optical 
data. Our data show weak 2.5 cr, 4.2 cr, and 2.1 cr features at the wavelengths of the Cr II 
A 2056, Cr II A 2062, and Zn II A 2062 lines. The strong feature appearing near 200 km 
s"^ in the Cr II A 2056 fine cannot be Cr II, because it is not seen in the Cr II A2062 fine 
(which is lower in / value by a factor of only 1.3) even though the SNR in the continuum 
is comparable. Furthermore, the velocity of this component would be too high compared to 
the redshifts of the Mg II and 21-cm absorption fines. We take this feature to be a Lyman-a 
forest line. 
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The only other metal lines available in the hterature for this DLA are the Mg II A A 2796, 
2803 lines at z = 0.2378 in lower resolution spectra. The 21-cm absorption line detected 
by Kanekar & Chengalur (2001) at z — 0.23780 ± 0.00002 shows a single component with a 
FWHM of 7.7 km s~^, indicating 621cm = 4.6 km s~^ for cold neutral hydrogen. 

We carried out profile fits using simultaneously all the lines of Zn II and Cr II that were 
included in our data, i.e., Zn II A 2062 and Cr II AA 2056, 2062, 2066. The dashed curves 
in Fig. 11 show the profiles for the best-fit N, b, and v parameters. The estimate of the b 
value obtained from the fits is 15.0 ± 7.7 km s^^. In Table 8, we list the derived Zn II and 
Cr II column densities, and the corresponding abundances. The formal estimate of the best 
fit Zn II column density in this DLA has a relatively large total error; therefore we list it 
only as a 3 cr upper limit. Also listed in Table 8 are results of single-component fits to the 
individual detected fines for several (fixed) choices of 6, i.e., 5, 10, 15, and 20 km s~^. 

For illustrative purposes we show, in Fig. 12, simulations of the Zn II A 2062 line for 
solar, 32% solar, and 10% solar metallicities ([Zn/H] = 0.0, -0.5, and -1.0), assuming the 
best- fit b and v values, b — 15.0 km s~^ and v — 8.7 km s~^. The parameters of the Cr II 
A 2062 fine are held fixed in these simulations. These simulations support the profile fitting 
results that the metallicity appears to be low. 

Our results are summarized in Table 9. The errorbars on the abundance estimates 
include photon noise and the uncertainties in b value, central velocity, and continuum place- 
ment. The [Zn/H] value for the DLA in Q0827-I-243 is not well constrained because of the 
very low SNR in the spectrum of this quasar. The [Zn/H] values appear to be low for the 
DLAs toward Q0738+313 and Q0952+179. 

4. The Metallicity-Redshift Relation for DLAs 

To examine the implications of our data for the mean metallicity-redshift relation of 
DLAs, we combined our data with other Zn measurements at higher redshifts. While doing 
this, we included only the "classical DLAs", i.e., systems with log Nm > 20.3. Furthermore, 
we excluded systems with radial velocities within 3000 km s~^ from the quasar emission 
redshifts, since these associated systems may not have similar abundances to those of the 
intervening systems at the same redshifts. The combined sample contains a total of 87 DLAs 
with Zn detections or limits, in the redshift range 0.09 < z < 3.90, and is summarized in 
Table 10. The sample includes our HST data, MMT data for 4 intermediate redshift DLAs 
recently obtained by us (Khare et al. 2004), and other data from the literature, as detailed in 
Table 10. While combining these measurements from different studies, we have normalized 
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them all to the same set of oscillator strengths and solar abundances that we have used for 
our HST data. Our HST observations have provided 3 new Zn measurements dX z < 0.39 
where no data existed before. This sample, together with our recent data at 0.6 < z < 1.5 
(Khare et al. 2004), has doubled the z < 1.5 sample and increased the z < 1 sample by a 
factor of 2.5. Most of the other recent observations have increased the samples at z >2. 

We now examine the evolution of the global Zn metallicity for these 87 DLAs. Of 
these 87 systems, 51 are detections and 36 are limits. One way to deal with the limits is to 
consider that the true Zn column density for a system with a limit is somewhere between 
zero and the upper limit value. With this in mind, we first constructed two extreme samples, 
a "maximum-limits" sample, where we treated the Zn limits as detections, and a "minimum- 
limits" sample where we treated the Zn limits as zeros. For any individual system, these 
cases cover the full range of possible values the Zn column densities can take. For the sake 
of pictorial illustration, we divided each of these samples into six redshift bins with 14 or 
15 absorbers each. Figure 13 shows the logarithm of the iV(H I)-weighted mean metallicity 
relative to the solar value for these six bins, plotted as a function of redshift. Each point 
is plotted at the median redshift in its bin. Horizontal bars for each bin extend from the 
lowest DLA redshift to the highest DLA redshift within that bin. The top panel shows 
the "maximum-limits" case. The filled circles in the bottom panel show the "minimum- 
limits" case with hmits treated as zeros. The vertical error bars in both panels show the 1 
(7 uncertainties in the A^ni-weighted mean metallicities, and include both the measurement 
uncertainty in the individual Zn II and H I column density measurements, and the sampling 
uncertainty. Further descriptions of the statistical procedures can be found in Kulkarni & 
Fall (2002). Clearly, Fig. 13 shows that the A^ni-weighted mean metaUicity of DLAs does 
not rise fast enough to reach solar or near-solar metallicities at 2; = 0. 

The above approach is perfectly adequate for low redshifts, which is the main redshift 
range of interest for our paper, since the number of Zn limits is relatively small at low 
redshifts. However, at high redshifts, there are a large number of Zn limits. At the same 
time, for most high- 2; systems with Zn limits, measurements of other elements are available. 
We therefore examined how much the results of the "minimum-limits" sample would change 
if we use information available from other elements for the systems with Zn limits. To do 
this, we created a third sample, using conservative constraints on the metallicities imphed by 
the information for the other elements (detections or lower limits) available in the literature 
(Lu et al. 1996; Ledoux et al. 1998; Ledoux, Bergeron, & Petitjean 2002; Prochaska & 
Wolfe 1999; Prochaska et al. 2001a; 2002; 2003a; 2003b; 2003c; and references therein). 
For 23 of the 36 systems with Zn upper limits, we estimated lower limits for Zn using 
measurements of the weakly depleted element Si, and assuming the total (gas + solid phase) 
relative abundance Si/Zn to be in solar proportion. For one system, we set the Zn abundance 
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to the abundance of the undepleted element S, assuming solar S/Zn ratio. For most of the 
remaining systems without Si, S, or O measurements, we set the constraints on Zn using 
measurements of Fe (8 systems) or Cr (2 systems). While setting the Zn constraints using 
Fe or Cr, we conservatively added 0.4 dex to correct for a typical amount of depletion in 
DLAs, unless other elements suggested no depletion. Finally, for two systems, no other 
measurements could be used to set a lower limit on Zn. In these cases we set the lower limits 
on Zn abundances to 0.01 solar values, since DLA metaUicities are seldom found to be below 
0.01 solar. The unfilled circles in the bottom panel of Fig. 13 show the mean metaUicity 
in each redshift bin for this modified minimum-limits sample with the information on other 
elements used in cases of Zn hmits. The two bins at z > 2.8 show some change as a result 
of using the other elements, but there is very little or no change in the other four bins at 
z < 2.8. Thus our conclusion about the low mean metallicity at low redshifts is not sensitive 
to whether or not information from other elements is used in cases of Zn limits. 

We now quantify the shape of the metallicity-redshift relation using several statistical 
tests outlined in Kulkarni & Fall (2002). Tabic 11 lists the results for the binned metallicity- 
redshift relation for the three samples considered above, in columns 4, 5, and 6. The log- 
arithm of the A^Hi-weighted mean metallicity in the lowest redshift bin 0.09 < z < 1.37 is 
—0.86 ± 0.11 for the maximum-limits sample, —1.01 ± 0.14 for the minimum-limits sample 
using Zn only (i.e. treating Zn limits as zeros), and —0.93 ± 0.10 for the minimum-limits 
sample using information from other elements in cases of Zn limits. For the 2; < 0.6 range 
not accessible to ground-based telescopes, the logarithm of the A^ni-weighted mean metallic- 
ity is -0.78 for the maximum-limits sample and -0.99 for the minimum-limits sample using 
information from other elements in cases of Zn limits. 

We next fit the data with a straight hne relation of the form log (Z/Zq) — log (Zo/Zq) + 
Bz, and hst the results for the parameters in Table 12. For the binned data (with six bins in 
the range 0.09 < z < 3.90), the linear regression slope B of the logarithm of the metaUicity 
vs. redshift relation is —0.18 ± 0.06 for the maximum-limits sample, —0.22 ± 0.08 for the 
minimum-limits sample with Zn only (i.e. with Zn limits treated as zeros), and —0.23 ±0.06 
for the minimum-limits sample using constraints from other elments in cases of Zn limits. The 
linear regression intercept of the metallicity-redshift relation (which corresponds to predicted 
present-day metallicity) is —0.74 ±0.15 for the maximum-limits sample, —0.75 ±0.18 for the 
minimum-limits sample with Zn only, and —0.71 ±0.13 for the minimum-limits sample using 
constraints from other elments in cases of Zn limits. Of course, as mentioned by Kulkarni & 
Fall (2002; see their footnote 3), unhke the means in individual bins, the slopes and intercepts 
thus determined do not necessarily bracket the true ranges of the slope and the intercept. 
This is because, in principle, all the true values in cases of upper limits could be closer to 
the limits at high redshifts, and closer to zero at low redshifts. Nevertheless, these cases are 
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indicative of the ranges the slope and the intercept are hkely to take. 

Another approximate, although not rigorously justified, way to get an indication of the 
slope and intercept of the metallicity-redshift relation from our maximum and minimum 
limits cases is to take the mean of the logarithmic metallicity estimates for the maximum 
limits and minimum limits cases in each redshift bin. Column (7) of Table 11 lists the means 
of the logarithmic metallicity estimates for the maximum limits case (column 4) and the 
minimum limits case with information from other elements in cases of Zn limits (column 6). 
The error bars on these mean values include the uncertainty in the values for the maximum 
and minimum limits cases and the standard deviation of these two values with respect to 
their mean value. The error bar on the mean value is large for the highest redshift bin 
(3.28 < z < 3.90), which contains mostly hmits. The hnear regression slope and intercept of 
the global metallicity-redshift relation obtained on using these mean logarithmic metallicity 
estimates in each bin are —0.20 ± 0.07 and —0.72 ± 0.16, respectively. 

We have also repeated several other calculations of the metaUicity-redshift relation ex- 
cluding the highest redshift bin (where there are mostly hmits rather than detections of 
Zn); using an unbinned A'"(HI)-weighted nonlinear approach; or using survival analysis, 
following the methodology of Kulkarni & Fall (2002). In the survival analysis calculations, 
we computed the Kaplan-Meier (K-M) estimator, using the astronomical survival analysis 
package ASURV rev. 1.2 (Isobe & Feigelson 1990) which implements the methods of Feigel- 
son & Nelson (1985). The error bars on the survival analysis estimates include the sampling 
and measurement uncertainties in the column densities. The results of all these calculations 
are listed in Tables 11 and 12. 

Finally, to assess the implications of our HST data alone for the statistics of the 
metallicity-redshift relation, we also repeated our analysis for the 83 systems left after tak- 
ing out the other new data used here for 4 intermediate-redshift systems from Khare et al. 
(2004). These calculations yielded essentially identical estimates of the metallicity-redshift 
relation. For example, for the maximum-limits sample, the A^Hi-weighted mean metallcity 
in the lowest redshift bin changes by only Ri 0.03 dex. The intercept and slope of the 
metaUicity-redshift relation also change by only < 0.03. 

The various tests described above give slightly different values, but they all confirm 
that the slope of the metallicity-redshift relation is small, but probably not zero, and the 
intercept is low. In other words, the global metaUicity of DLAs does not seem to rise to the 
solar value at low redshifts. This just confirms the visual impression from Fig. 13. 

Our estimates of the metallicity-redshift relation are consistent with previous results 
of Kulkarni & Fall (2002) and Prochaska et al. (2003b). We have not included metaUicity 
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estimates based on X-ray absorption in our analysis here because the X-ray data are expected 
to reflect mainly the abundance of O, and there could be a systematic difference between 
[Zn/H] and the [X/H] derived from the X-ray absorption data. ^ 

To summarize, our analysis suggests that the global metaUicity of DLAs stays low at 
low redshifts, and evolves slowly at most. We stress, however, that the current samples are 
still small, and further observations, particularly at the low and intermediate redshifts are 
necessary to decrease the uncertainties in the shape of the metallicity-redshift relation. 

We now compare our results with some models of cosmic chemical evolution. The 
short-dashed, long-dashed, solid, and short-dash-long-dashed curves in Figure 13 show, re- 
spectively, (a) the mean interstellar metallicity in the closed-box and outflow models of Pei 
& Fall (1995); (b) the mean Zn metallicity in the Malaney & Chaboyer (1996) model with a 
constant Zn yield, a slope of —1.70 for the initial mass function (IMF), and the evolution of 
the neutral gas density as computed by Pei & Fall (1995); (c) the mean interstellar metallic- 
ity in the Pei et al. (1999) model with the optimum fit for the cosmic infrared background 
intensity; and (d) the mean metallicity of cold interstellar gas in the semi-analytic model of 
Somerville et al. (2001). Models (a) and (c) agree with the data at high redshifts, but evolve 
too rapidly toward solar metaUicity to be consistent with the data at low redshifts. Models 
(b) and (d) have lower slopes but predict too high a metallicity at all redshifts. None of the 
models agrees completely with all the observations. 



5. Discussion 

The main new result from our study is that, in all the 3 DLAs for which our observations 
offer meaningful constraints on [Zn/H], the metallicities appear to be substantially sub-solar, 
and that there is at most only a weak evolution in the global mean metallicity with redshift. 

To understand this puzzling result, we now explore (a) whether our results are affected in 
any way by the use of Zn as a metallicity indicator, and (b) how our results fit in with 



^Prochaska ct al. (2003b) reported a slope of —0.25 ± 0.07 for the metallicity-redshift relation in the 
redshift range 0.5 < z < 4.7., a detection of evolution in the global mean metallicity at > Str level. Their 
analysis included the X-ray measurement from Junkkarinen et al. (2004) for the z = 0.52 DLA toward 
Q0235+164. This system, because of its low redshift and very high H I column density (log TVhi = 21.80), 
has a very strong influence on the mean metallicity in the lowest redshift bin, and therefore on the overall 
slope of the mean metallicity-redshift relation. When Kulkarni & Fall (2002) included X-ray absorbers in their 
analysis, they also found evolution at > Str level (see their foonote 4). For most of their analysis, Kulkarni & 
Fall (2002) left out X-ray absorbers since the X-ray measurements primarily trace the O abundance, which 
could differ systematically from the Zn abundance. 
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the morpohological information available for these DLAs. Finally, we discuss the role of 
potential selection effects caused by dust obscuration in the DLAs. 

5.1. Zinc as the Metallicity Indicator 

Could the apparent low metaUicity for most of our target DLAs be an artifact of the use 
of Zn as the metallicity indicator? To explore this question, we now examine the advantages 
and disadvantages of using Zn. 

The conclusion that [Zn/Fe] in Galactic halo and disk stars for metaUicities [Fe/H] 
> —3 that followed from earher work (Sneden, Gratton, & Crocker 1991) has been strongly 
supported by recent studies (e.g., Mishenina et al. 2002; Cayrel et al. 2004; Nissen et 
al. 2004). There is some evidence for a small overabundance of Zn by ~ 0.1 — 0.2 dex in 

metal-poor disk stars and halo stars with [Fe/H] < —2 (e.g., Prochaska et al. 2000; Bihain 
et al. 2004; Cayrel et al. 2004). Chen ct al. (2004) have also recently reported [Zn/Fe] 
~ 0.1 — 0.15 dex in disk stars with [Fc/H] > —0.9. However, in all of these studies, there 
are systematic uncertainties in [Zn/Fc] at the level of ^ ±0.1 dex, and overall Zn tracks 
Fe quite closely for [Fe/H] > —2. These observations are not completely understood in 
nucleosynthetic models of massive stars. The models are complex and the value of [Zn/Fe] 
depends on several parameters such as the position of the mass cut, the explosion energy, 
the degree of neutronization, and the consequences of a neutrino-driven wind (Woosley & 
Weaver 1995; Chicffi & Limongi 2002; Umeda & Nomoto 2002; Woosley, Heger, & Weaver 
2002; Fenner, Prochaska, & Gibson 2004). Despite these uncertainties in the theoretical 
understanding, the Galactic stellar data support the use of Zn as a tracer of metallicity. 

There are some observational difficulties in measuring the Zn abundances in DLAs. The 
Zn II AA2026, 2062 lines are weak and sometimes hard to detect. On the other hand, for 
that same reason, they are usually unsaturated. At 2; > 3, the Zn II lines lie in spectral 
regions affected by night sky lines, and thus become difficult to detect. At the same time, 
the long wavelengths of the Zn II lines also make them accessible in near-UV and optical 
spectra for the range z <3 which spans ~ 90% of the cosmic history. 

While some studies have suggested that ionization corrections for Zn could be significant 
(e.g., Howk & Sembach 1999), others have given convincing arguments that such corrections 
are negligible. For example, Vladilo et al. (2001) have pointed out that the ionization 
corrections are uncertain because of uncertainties in the Zn recombination coefficients and 
ionization cross sections. They have provided evidence that the Zn ionization corrections 
could be lower, and concluded that DLA metallicity estimates based on [Zn/H] are unlikely 



-20- 



to be significantly affected by ionization effects. Indeed, at the large H I column densities 
(log A^Hi = 20.9 — 21.3) of the 3 DLAs in our sample where we find low [Zn/H], the ionization 
correction for [Zn/H] is expected to be at most about 0.3 dex, and probably smaller for more 
accurate atomic parameters (see, e. g.. Fig. 16 of Vladilo et al. 2001). 

Thus, while Zn, hke all other elements, has both advantages and disadvantages, it is still 
an appropriate choice as a metallicity indicator in DLAs (sec, e.g., Bihain et al. 2004; Nissen 
et al. 2004). Furthermore, the various uncertainties mentioned above are not large enough 
to make significant difference to our results regarding a low global mean metallicity -0.7 
to -1.0 dex) at low rcdshifts. Overall, we conclude that our main result is not affected much 
by the use of Zn as the primary metallicity indicator. 

5.2. Compcirison with Morphological Information 

Could the low metallicities in our low-redshift DLAs be understood in terms of the 
nature of the underlying galaxies? In fact, the low metallicities are consistent with the 
morphological information available on the candidate galaxies believed to give rise to these 
DLAs. 

The z = 0.0912 DLA toward Q0738-I-313 has been tentatively identified as associated 
with fuzz surrounding the quasar, and classified as a low surface brightness galaxy. This 
galaxy has an estimated luminosity of 0.08L^ and impact parameter < 3.6 kpc for 
rim ^ l,flA = 0, and Ho-^65h km s"^ Mpc^^ (Cohen 2001; Turnshek et al. 2001). The 
z — 0.2212 DLA toward the same quasar has been spectroscopically identified as a dwarf 
spiral galaxy, with B and K-band luminosities of O.IL^ and O.IL^, respectively, and impact 
parameter of 20 h'^ kpc (Cohen 2001; Turnshek et al. 2001). 

From ground-based images in the B, R, J, and K bands, Rao et al. (2003) have suggested 
that the best candidates for the z — 0.2378 DLA toward Q0952-I-179 are two dwarf low 
surface brightness galaxies with a combined luminosity of 0.02L^, and impact parameter < 
4.5 kpc. However, these candidate objects have not yet been confirmed spectroscopically. 
In any case, the DLA in this field cannot be an L* galaxy, as none of the luminous galaxies 
in the field have the same redshift as the DLA (Bergeron & Boisse 1991). 

For Q0827-I-243 {zabs — 0.5247), where we have only a loose upper limit on the metal- 
licity, a galaxy at 2; = 0.5258 has been spectroscopically confirmed at an angular separation 
of 6", corresponding to a projected separation of 34 kpc. This galaxy has B and K-band 
luminosities of 0.8L^ and 1.2L^, respectively. Rao et al. (2003) have suggested a "disturbed 
spiral" morphology for this galaxy. Unfortunately, our Zn data do not yield a firm deter- 
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mination of the metallicity in this DLA. The Fe abundance is low ([Fe/H] = —1.02 ± 0.05; 

Khare et al. 2004); but this could be partly because of depletion of Fe on dust grains, and 
the total (gas + solid phase) metallicity could be higher. In any case, at the large impact pa- 
rameter for this absorber, the metallicity is likely to be low, for a typical abundance gradient 
of ~ 0.05 — 0.1 dex kpc^^ observed in the Milky Way (e.g., Friel & Janes 1993). 

5.3. Potential Selection Effects 

Our observations of low metallicities in low-redshift DLAs are surprising because the 
mean mass-weighted interstellar metallicity of nearby galaxies is found to be close to solar 
(e.g., Kulkarni & Fall 2002). Thus, samples derived from absorption-line studies of optically 
selected quasars seem to be biased toward less metal-rich galaxies than would be expected 
from an H I emission-selected sample of galaxies. 

The low metallicities of our low- 2; DLAs also appear to contradict the expectations of 
a solar or near-solar metallicity at z = from most cosmic chemical evolution models based 
on the global star formation history of galaxies. Our observations thus suggest that the 
current DLA samples, especially those at low redshifts, might not be representative of the 
global galaxy population. If the observed DLAs have low metallicities at all redshifts, it is 
possible that the DLA galaxies are dominated by dwarf or low surface brightness galaxies. 
Some earher studies have suggested that DLAs could arise in dwarf galaxies (e.g., York et 
al. 1986, Matteucci et al. 1997). Indeed, some imaging studies have suggested DLAs at 
low redshifts to be predominantly dwarf or low surface brightness galaxies (e.g., Rao et al. 
2003). In fact, as discussed in section 5.2, three out of the four DLAs in our sample are 
tentatively identified as dwarf or low surface brightness galaxies (although only one of them 
is spectroscopically confirmed). Again, taken at face value, this is very surprising because 
the H I emission studies of galaxies at z = show that most of the H I at the present 
epoch resides in large spirals (e.g., Rao & Briggs 1993; Zwaan, Briggs, & Sprayberry 2001). 
Of course, the identification of low- 2; DLAs with dwarf or low surface brightness galaxies 
is not always certain, since many imaging studies (e.g., LeBrun et al. 1997; Rao et al. 
2003) have not always confirmed the redshifts of the candidate galaxies. Chen & Lanzctta 
(2003) have obtained more extensive photometric redshift data for DLA fields and suggested 
that the contribution of dwarf galaxies to DLAs at 2; < 1 is not significant. Rosenberg & 
Schneider (2003) have also suggested that the luminosity function of low-redshift DLAs may 
be fairly fiat over several orders of magnitude. Thus, there seems to be no full consensus 
about the nature of low-redshift DLAs. Clearly, it is necessary to obtain further imaging 
studies and redshift identifications of low-redshift DLAs and take into account the luminosity 
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uncertainties of the DLA galaxies. 

Our observations of low metallicities in low-redshift DLAs suggest that these galaxies 
have had less nucleosynthetic processing and hence low time-averaged star formation rates 
(SFRs). On the other hand, observations of the UV luminosity density of galaxies (e.g. Lilly 
et al. 1996; Madau et al. 1996) indicate that the global average SFR in galaxies was high 
at z > 1.5. In fact, it is this high global SFR at high redshifts that drives the metallicity to 
rise up to solar or near-solar levels by the present epoch in most cosmic chemical evolution 
models. There arc also indications from C II* absorption that some DLAs may have high 
SFRs and could be efficient producers of heavy elements (e.g., Wolfe ct al. 2003). If the 
observed low-redshift DLAs are largely metal-poor, then where are the more metal-rich low- 
z absorbers expected to have near-solar metaUicities? Why are the absorption studies not 
detecting more luminous and metal-rich DLAs even at low redshifts? Such a selection effect 
could be caused by dust obscuration in the DLAs. This is because the DLA galaxies with the 
highest column densities of metals arc likely to be those with the highest column densities 
of dust, and these may obscure background quasars to such a degree that some of them are 
omitted from optically selected samples (Fall & Pei 1993; Boisse et al. 1998). As a result 
of obscuration, the observed mean metallicity in the DLA galaxies, and hence the data 
points in Figure 13, could lie systematically below the true mean metallicity. The model 
predictions plotted in Fig. 13, however, refer to the "true" mean metallicity (not corrected 
for dust obscuration). Indeed, if the true mean metallicity and dust content increase with 
decreasing redshift, more DLAs could be missing at lower redshifts. This would flatten the 
metallicity-redshift relation, as possibly shown by our observations. 

In principle, a comparison of metal abundances in DLAs from optical vs. radio-selected 
quasars can help to quantify the dust selection effect. At high redshifts, EUison et al. (2001) 
found slightly more DLAs in the foreground of radio-selected and optically faint quasars than 
in optically selected and optically bright quasars. This result is sometimes misinterpreted as 
evidence against dust selection effects. However, the radio-selected samples arc still small, 
and the uncertainties correspondingly large. The difference between the true (unobscured) 
and observed column density distributions /(iVni) of DLAs is expected to be only < 50% for 
log A^Hi < 21 (see Fig. 1 of Pei & Fall 1995). The Ellison et al. (2001) points, within the error 
bars, are in fact consistent with the predictions of Fall & Pei (1993), Pei & Fall (1995), and 
Pei et al. (1999). Recent results of EUison et al. (2004) for Mg II systems with 0.6 < ^ < 1.7 
also appear to be consistent, within the large error bars, with the predictions of Fall & Pei 
(1993), Pei & Fall (1995), and Pei et al. (1999), although exact comparisons of /(A^m) of 
DLAs can only be made after the DLA nature of the Mg II systems is confirmed with UV 
spectra. In a study of reddening in quasars with and without DLAs, Murphy & Liske (2004) 
have recently suggested that the extent of dust reddening caused by high-redshift DLAs may 
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be small at z ~ 3, but the effects at low redshift are likely to be larger. Indeed, Vladilo 
(2004) has reported that the dust-to-metals ratio in DLAs, as inferred from the depletion of 
Fe, increases with decreasing redshift. Abundance studies for a very large sample of DLAs 
in radio-selected and optically faint quasars would help to improve the constraints on the 
extent of the dust obscuration bias. 



6. Summciry 

We have obtained HSTSTIS UV spectra of three quasars with four DLAs in the redshift 
range 0.09 < z < 0.52, to constrain the abundances of Zn, Cr, and Fe in these absorbers. 
These observations have provided the first constraints on Zn abundances in DLAs with 
z < 0.4. In all the 3 DLAs for which our observations offer meaningful constraints on the 
Zn abundance, the metallicities appear to be substantially sub-solar (< 10 — 20% solar). 
We have combined our results with higher redshift data from the literature to estimate the 
global mean metallicity-redshift relation for DLAs. We find that the global mean metaUicity 
of DLAs shows at most a slow increase with decreasing redshift and does not appear to rise 
up to solar or near-solar values at ^ = 0. On the other hand, most cosmic chemical evolution 
models, based on the global star formation history inferred from the emission studies of 
galaxies, predict the global interstellar metallicity to be solar or near-solar at z = 0. This 
suggests that the DLAs, especially those at low redshifts, do not trace the general population 
of galaxies represented by the global star formation history of the universe. The low A^(H 1)- 
weighted mean metallicity of low-redshift absorbers also appears to contradict the fact that 
the mass-weighted mean interstellar metallicity of nearby galaxies is observed to be close to 
solar. 

This weak evolution in the DLA global mean metallicity could be explained by the 
fact that our sample seems to be dominated by dwarf or low surface brightness galaxies. It 
is possible that current DLA samples, especially those at low redshifts, are biased against 
more enriched galaxies because the latter may cause more dust obscuration of the background 
quasars. We note, however, that the low-redshift samples are still small, and the fact that 
most of the DLAs in our sample are metal-poor could be an effect of small number statistics. 
Metallicity measurements of more low-redshift DLAs are necessary to further improve the 
statistics of the mean metallicity-redshift relation and to study the scatter around the mean. 
Observations of low-redshift DLAs toward optically faint quasars arc especially necessary in 
the future to improve the constraints on the mean metallicity-redshift relation for DLAs, 
and to quantify potential dust selection effects. 
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Fig. 6.— Velocity plot of the Fe II, Zn II, and Cr II data for the z = 0.0912 DLA toward 
Q0738+313. The zero point of the velocity scale refers to the absorption redshift z = 0.0912 
inferred from the 21-cm and optical data. The dashed curves show the profiles for the best-fit 
A^, b, and v parameters obtained by fitting multiple lines simultaneously. In the case of Zn 
II, the profiles correspond to the 3 a upper limit adopted. 
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Fig. 7.— Comparison of Zn II A2026 data for the z = 0.0912 DLA toward Q0738+313 with 
simulated Voigt profiles for [Zn/H] = -1.0, -0.5, and 0.0, shown by solid (red), short-dashed 
(green), and short-dash- long-dashed (blue) curves, respectively. All three curves assume the 
best-fit b and v values obtained by fitting multiple lines simultaneously. 
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Fig. 8.— Velocity plot of the Zn II and Cr II data for the z = 0.2212 DLA toward Q0738+313. 
The zero point of the velocity scale refers to the absorption redshift z = 0.2212 inferred from 
the 21-cm and optical data. The dashed curves show the profiles for the best-fit A^, b, and v 
parameters obtained by fitting multiple lines simultaneously. In the case of Zn II, the profiles 
correspond to the 3 a upper limit adopted. 
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Fig. 9.— Comparison of Zn II A2062.664 data for the z = 0.2212 DLA toward Q0738+313 
with simulated Voigt profiles for [Zn/H] = -1.0, -0.5, and 0.0, shown by sohd (red), short- 
dashed (green), and short-dash-Iong-dashed (blue) curves, respectively. All three curves 
assume the best-fit values of b, v, and Cr II column density obtained by fitting multiple lines 
simultaneously. 
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Fig. 10.— Velocity plot of the Zn II and Cr II data for the z = 0.5247 DLA toward 
Q0827+243. The zero point of the velocity scale refers to the absorption redshift z = 0.5247 
inferred from the 21-cm and optical data. The dashed curves show profiles corresponding to 
the 3 cr upper limits adopted for 6 = 10 km s~^. 
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Q0952+179 zabs = 0.2378 




Fig. 11. — Velocity plot of the Zn II and Cr II data for the z = 0.2378 DLA toward 
Q0952+179. The zero point of the velocity scale refers to the absorption redshift z = 0.2378 
inferred from the 21-cm and optical data. The dashed curves show the profiles for the best-fit 
A^, b, and v parameters obtained by fitting multiple lines simultaneously. In the case of Zn 
II, the profiles correspond to the 3 a upper limit adopted. 
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Fig. 12.— Comparison of Zn II A2062 data for the z = 0.2378 DLA toward Q0952+179 
with simulated Voigt profiles for [Zn/H] = -1.0, -0.5, and 0.0, shown by sohd (red), short- 
dashed (green), and short-dash-Iong-dashed curves (blue), respectively. All three curves 
assume the best-fit values of b, v, and Cr II column density obtained by fitting multiple lines 
simultaneously. 
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Fig. 13. — The global metallicity-redshift relation deduced from damped Ly-a absorbers. 
The circles show the logarithm of the A^(H I)-weighted mean Zn metallicity relative to the 
solar value vs. redshift for the sample of 87 DLAs based on our HST observations and 
the literature. The filled circles in the top (a) and bottom (b) panels show, respectively, 
the A^(H I)-weighted mean metallicity for the maximum-limits and minimum-limits cases 
using Zn only. The unfilled circles in the bottom panel show the A^(H I)-weighted mean 
metallicity for the modified minimum- limits sample, using other elements to constrain the 
metallicities in cases of Zn limits. Vertical error bars on the filled circles in both panels 
denote 1 a uncertainties in the logarithm of the A^(H I)-weighted mean metallicity. Data 
points are plotted at the median redshift in each bin. Horizontal bars for each bin extend 
from the lowest DLA redshift to the highest DLA redshift within that bin. The short-dashed 
(blue), long-dashed (green), solid (red), and short-dash-long-dashed (magenta) curves show, 
respectively, the "true" mean metallicity (not corrected for dust obscuration) expected in 
the cosmic chemical evolution models of Pei & Fall (1995), Malaney & Chaboyer (1996), Pei 
et al. (1999), and Somerville et al. (2001). 
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TABLE 1 
Target DLAs 



QSO 




^abs 


log Nm 


Rcf.t 


Q0738+313 
Q0738+313 
Q0827+243 
Q0952+179 


0.635 
0.635 
0.939 
1.472 


0.0912 
0.2212 
0.5247 
0.2378 


21.181°:°^ 

20.301°:°^ 

91 09+0.05 


1,2 
1,3 
4,5 
4,5 



I References: 1. Rao & Turnshck (1998); 2. Chengalur & Kanekar (1999); 3. Lane et al. 
(1998); 4. Rao & Turnshek (2000); 5. Kanekar & Chengalur (2001). 
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TABLE 2 
Details of HST STIS Observations 



QSO 


Detector 


Grating 


R 


^ccnt 


^cxp ip) 


Q0738+313 


NUV-MAMA 


G230M 


12,100 


2176 


13,388 


Q0738+313 


NUV-MAMA 


G230M 


12,500 


2257 


11,224 


Q0738+313 


NUV-MAMA 


G230M 


13,900 


2499 


13,613 


Q0827+243 


CCD 


G230MB 


10,400 


3115 


12,232 


Q0952+179 


NUV-MAMA 


G230M 


14,300 


2579 


24,346 



-45 - 



TABLE 3 
Line Lists 



QSO 


No. 






p 


ID 








(A) 


(A) 


(A) 






Q0738+313 


It 


2155.47 


0.140 


0.039 








2t 


2187.08 


0.157 


0.043 








3 


2211.03 


0.092 


0.027 


Zn II 2026 


0.0912 




4t 


2216.58 


0.157 


0.039 








5 


2243.99 


0.096 


0.029 


Cr II 2056 


0.0912 




6 


2260.75 


0.180 


0.033 


Fe II 2260 


0.0000 




7 


2467.07 


0.109 


0.028 


Fe II 2260 


0.0912 




8 


2504.76 


0.403 


0.036 


C IV 1548 


0.6179 




9 


2508.86 


0.234 


0.030 


C IV 1551 


0.6178 




10 


2518.49 


0.061 


0.020 


Cr II 2062 


0.2212 




11 


2519.13 


0.066 


0.017 


Zn II 2062 


0.2213 


Q0952+179 


1 


2545.32 


0.070 


0.020 


Cr II 2056 


0.2378 




2 


2546.78 


0.234 


0.032 


Ly-a 


1.0950 




3 


2552.69 


0.100 


0.024 


Cr II 2062 


0.2378 




4 


2553.18 


0.052 


0.021 


Zn II 2062 


0.2378 




5 


2560.01 


0.060 


0.019 


Ly-a 


1.1058 




6 


2564.45 


0.616 


0.037 


Ly-a 


1.1095 




7 


2566.22 


0.770 


0.033 


Ly-a 


1.1109 




8 


2571.61 


0.081 


0.025 


Ly-a 


1.1154 




9 


2576.78 


0.282 


0.036 


Mn II 2577 


0.0003 




10 


2586.60 


0.843 


0.042 


Fe II 2586 


0.0003 




11 


2590.69 


0.447 


0.044 


Ly-a 


1.1311 




12 


2594.27 


0.262 


0.039 


Mn II 2594 


0.0002 




13 


2600.21 


1.054 


0.043 


Fell 2600 


0.0003 




14 


2606.26 


0.218 


0.040 


Mn II 2606 


0.0002 




15 


2612.30 


0.198 


0.029 


Ly-a 


1.1489 




16 


2615.42 


1.224 


0.057 


Ly-a 


1.1514 



t No clear identification, but possibly C IV A 1548 lines at Zahs — 0.392, 0.413, and 
0.432 with the corresponding C IV A1551 hnes lost in noise. 



-46- 



TABLE 4 

Further Details on DLA Metal-line Features 



QSO 


^obs 




p 






ID 






^abs 




(A) 


(A) 


(A) 


(A) 


(A) 










Q0738+313 


2211.03 


0.092 


0.027 


0.019 


0.033 


Zn 


II 


2026 


0.0912 




2243.99 


0.096 


0.029 


0.020 


0.036 


Cr 


II 


2056 


0.0912 




2467.07 


0.109 


0.028 


0.018 


0.033 


Fe 


II 


2260 


0.0912 


Q0738+313 


2518.49 


0.061 


0.020 


0.009 


0.022 


Cr 


II 


2062 


0.2212 




2519.13 


0.066 


0.017 


0.008 


0.019 


Zn 


II 


2062 


0.2213 


Q0952+179 


2545.32 


0.070 


0.020 


0.019 


0.028 


Cr 


II 


2056 


0.2378 




2552.69 


0.100 


0.021 


0.012 


0.024 


Cr 


II 


2062 


0.2378 




2553.18 


0.052 


0.021 


0.013 


0.025 


Zn 


II 


2062 


0.2378 
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TABLE 5 

Column Densities and Abundances for the z = 0.0912 DLA toward Q0738+313 



Fitting 
Procedure 


h 

(km 


log iVcrll 

(cm~^) 


log iVznII 

(cm-2) 


log A^Fell 

(cm~^) 


[Cr/H] 


[Zn/H] 


[Fe/H] 


N ^^v varied^; 
multiple lines 


11.2 
±6.7 


1 'X 08+0-l^ 


< 12.66 


15.021°:- 


1.00_o.23 


< -1.14 


-1-631^:}^ 


h fixed; 
single lines 


5.0 
10.0 
15.0 
20.0 


1 3 56+°-2'^ 

iO.OU„Q 5g 

13.491°:^^ 
13.48l°il 
13.501°;}^ 


13.061°:! 
12.841°:!^ 

12.811°:- 
12.811°:- 


15.221°:^! 

15.061°:^;^ 

15.061°:- 
15.061°:- 


1 97+0.25 
-'-•^'-0.56 
_1 34+017 

1 qr+0.15 

1 OO+0.17 
l.OO_o,20 


^- ' '-'-0.49 

-0 97+°-i'^ 

^•^ ' -0.20 

-1 00+°-^''^ 
-1 00+°-^^ 

-'-•'""-'-0.18 


-'-•^'J-0.41 

i.oy_Q 21 

-'-•'-'^-0.18 
-•-•3^-0.17 



f For the fits with iV, 6, f varied, the error bars on column densities and abundances 
include uncertainties from Poisson statistics and uncertainties in h value, central velocity, 
and continuum placement. The limits are 3 cr upper limits. 



-48- 



TABLE 6 

Column Densities and Abundances for the z = 0.2212 DLA toward Q0738+313 



Fitting Procedure 


b 

(km s^-*^) 


log Ncrll 

(cm^^) 


log iVznII 

(cm~^) 


[Cr/H] 


[Zn/H] 


N,b,v varied^; multiple lines^^ 


12.2 ± 7.5 


10.ll_o.24 






-0.70l°:i^ 


b fixed; single lines 


5.0 
10.0 
15.0 
20.0 


1 q 99+0.16 
10.ZZ_o.25 

-l-O.ZU_o 20 
-|o 21+0.13 
-'-"^■■^-'--0.19 

13.23l°:i^ 


1 9 7c;+0.15 
IZ. ( O_o.24 

12.73l°:J3 
12.76l^:i? 
12.801°:}^ 


1 qq+0.18 
1.00_o 26 
1 qc,+0.16 
l.OO_02i 

_1 Q4+0.15 
-'-■'-•^-0.20 
1 q9+0.15 
l.OZ_o 20 


-0.801°:}^ 

-0.73l°:i^ 



f: For the fits with A^, b, v varied, the error bars on column densities and abundances 
include uncertainties from Poisson statistics and uncertainties in b value, central velocity, 
and continuum placement, ff: The noisy data at the positions of the A2026 and A2056 
lines are not included in these fits, but they give consistent results of log iVznu < 12.81 and 
log A^crii < 13.37, respectively for the best-fit b value of 12.2 km s~^. Nevertheless, to be 
conservative, we finally adopt an upper limit [Zn/H]< —0.70 for this system. 
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TABLE 7 

Column Densities and Abundances for the z = 0.5247 DLA toward Q0827+243 



Fitting Procedure 


h 

(km s^"*^) 


log A^Crll 

(cm^^) 


log iVznII 

(cm^^) 


[Cr/H] 


[Zn/H] 


h fixed; single lines^ 


5.0 


< 13.79 


< 13.23 


< -0.16 


< 0.30 




10.0 


< 13.51 


< 12.89 


< -0.44 


< -0.04 




15.0 


< 13.45 


< 12.83 


< -0.50 


< -0.10 




20.0 


< 13. 12 


< 12.80 


< -0.53 


< -0.13 



f The limits are 3 a upper limits. 



-50- 



TABLE 8 

Column Densities and Abundances for the z = 0.2378 DLA toward Q0952+179 



Fitting Procedure 


h 

(km s^"^) 


log A^Crll 

(cm^^) 


log A^Znll 

(cm^^) 


[Cr/H] 


[Zn/H] 


A^,6,f varied"''; 
multiple lines 


15.0 ±7.7 


13.32l°:}2 


< 12.93 


-1.651°;}^ 


< -1.02 


h fixed; single lines 


5.0 


13.30 ±0.20 


12.62 ±0.20 


-1.67 ±0.21 


-1.33 ±0.25 




10.0 


13.24 ±0.20 


12.62 ±0.20 


-1.73 ±0.21 


-1.33 ±0.25 




15.0 


13.25 ±0.20 


12.66 ±0.20 


-1.72 ±0.21 


-1.29 ±0.25 




20.0 


13.26 ±0.20 


12.71 ±0.20 


-1.71 ±0.21 


-1.24 ±0.25 



f : For the fits with iV, 6, v varied, the error bars on column densities and abundances 
include uncertainties from Poisson statistics and uncertainties in h value, central velocity, 
and continuum placement. The limits are 3(T upper limits. 
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TABLE 9 
Summciry of Abundances 



QSO 


^ahs 


[Zn/H]* 


[Cr/H]* 


[Fe/H] 


Ref.t 


Q0738+313 
Q0738+313 
Q0827+243 
Q0952+179 


0.0912 
0.2212 
0.5247 
0.2378 


< -1.14 

< -0.70 

< -0.04 

< -1.02 


1 rr+0.15 
i.OO_Q 23 

_1 44+0-17 

< -0.44 
-1 65+°-^'' 


-i-63lS:i^ 

-1.02 ±0.05 


1 
1 

1,2 
1 



*: Values obtained from fits with A^, 6, v varied simultaneously for multiple lines, 
t References: 1. This Work; 2. Khare et al. (2004). 
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TABLE 10 
Summciry of Zn Measurements in DLAs 



QSO 


Ref.^ 


^dbs 


fZn/Hl 

L / J 


iVHT/1020(cm-2) 


iV7nTT/10^2(cm-2) 




1 

1. 


0Q1 9 


< —1 14 

\ X . x^ 


1 c, fin -(- 9 nn 


< 4 




1 

± 


n 9919 




7 Qfl -L 1 An 

1 .yu zn X .^u 


< RQ 

U.Ut/ 




± 


n 9'^78 


< —1 09 


91 fin -L 9 Kfi 
zx.uu zn z.uu 


^ O.UU 


1 99Q-n91 


9 


U.Oc/UU 


-0 4^ 


CC «9 -L fl Q1 

tj.uz zn u.yx 


Q CI _|_ 9 00 

O.UX ZIZ Zi.OO 




1 
1 


^947 


^ 04 


9 flfi _|_ fi 9fi 
z.uu in u.zu 


< 7 7^ 




Q 
O 


fi81 Q 


< _1 9Q 
\ — X . zy 


9 09 _L fl QQ 

z.oz in u.yy 


_\ U.U X 




Zl 


fiQ99 


1 90 
X . zu 


1 Q c;f) _|_ 9 9 c; 
xy.ou zn z.Ztj 


0.^ t ZIZ X 1 






7404 


fiO 


1 n 00 + 7^ 
xu.uu m u. ( 


1 n 79 _L 1 9/1 

xu . 1 Z( ZIZ X . Zt^ 


OzL^zL-i-O'^q 


7 




QQ 
u . y y 


A Qfl -1- fi8 
^.yu zn u.uo 


9 14 -h "^9 

Zi.X^ ZIZ U.OZi 


1 797-1- f;'?n9 


c; 
u 


n Q44Q 


-0 ^9 


14 cfi -1- 4 oc 
x^.uu zn ^.ou 


1Q «9 -L 1 Q1 
XO.UZ ZIZ X.t^X 


n'?n9-99'^ 


7 


1 nriQ'^ 

X.UUc/O 


-0 ^4 


9 90 -1- fl c;q 
z.zy zn u.uy 


9 09 -L fl 14 

Zi.OZi ZIZ u.x^ 






1 n'^19 


-X .oy 


9« fin + 7 sn 

zu.uu in / .ou 


/ ^ u.^u 




Q 
o 


1 1 ^nn 


u.yu 


9 09 -L 1 fifi 


1 '^fi -1- 1 

X.OU ZIZ u.xo 




s 

o 








xu.x^ ZIZ u.uu 




Q 


X.O ( 


78 


Zr.ou HZ u.ou 


X. ( U ZIZ u.^x 




1 


1 4900 


1 ^8 
X . oo 


O^.U 1 ZIZ ^.ox 


OQ 1 n q7 
O.Ocf ZL \J.cf 1 






1 47Q0 


1 ^8 
- X . oo 


AO nn -\- R 00 

^Zr.UU ZIZ O.UU 


^.UO Zl X ( 




1 1 


X . uu x^ 


QQ 
-u.yy 


i .UO ZIZ u. xu 


Q -1- 07 

O.UO Zl U.U ( 




1 9 

xz 


1 77'^fi 


-1 ^"^ 
X .uo 


90 fin _|_ c fifi 

ZiU.uu ZIZ u.uu 


9 en -1- ^0 
z.uu ZL u.uu 




1 


1 77fi4 


-1 9fi 
X .zu 


1 c no -1- 1 49 

xu.uu ZIZ X.^Zi 


_|_ n 90 

0.'4:0 ZL U.ZO 




14 


1 8'i7'^ 


-1 ni 

-X .ux 


7 _|_ 1 

/ .c/^ ZIZ X .ou 


001 _|_ n 77 

O.OX ZL u. / / 


99'?0-J-09^ 


1 ^ 


1 8fi49 


n fi8 

L/.UO 


7 no _L 1 00 
1 .uo ZIZ X .00 


a OA -u n 41 

U.OX ZL U.TiX 




1 


1 8Q90 


8fi 
-u.ou 


Q Qo 4. n Q9 
O.c/O ZIZ U.c/Zr 


ooAinAa 
Z.O^ ZL U . XU 


990fi 1 QQ 


1 


1 Q90^ 


'?7 


^.UU ZIZ U. ( ^ 


90 -u n 1 7 

O.ZU ZL U. X ( 


1 1 ^7-1-01 4 


1 (S 


1 Q440 


1 '^4 

- X . 


UO.XU ZIZ x^.uu 


1 9 QO 9 9« 

XZ.OU ZL Z.ZO 


O^^l '?(S(S 


1 7 


1 q(S90 


on 

-U . X X 


0. xu ZIZ U .Uc/ 


10 47 -k 1 91 

XU.t:/ zl X.ZX 






x.y / oo 


n 79 


70 -1- 70 
U. / U ZIZ u. / u 


0.^1 Zl U.uu 


1 8f;n-i-44 


1 Q 


1 Qonn 

x.yyuu 


-U.Oo 


9 1 9 _|_ Q/1 

ZU.XZ ZIZ U.O'4: 


99 on _|_ c 90 
ZiZi.oa Zl u.zu 


1215+333 


13,15 


1.9990 


-1.25 


8 91 ± 1 38 


2.14 + 0.24 


0010-002 


8 


2.0250 


-1.16 


6.31 ± 1.47 


1.86 + 0.28 


0458-020 


13,15 


2.0395 


-1.15 


44.67 + 9.32 


13.61 + 0.63 


2231-00 


15 


2.0660 


-0.73 


3.63 + 0.84 


2.90 + 0.15 


0049-283 


20 


2.0713 


< -1.04 


3.10 + 0.40 


< 1.22 


2206-199 


13 


2.0762 


< -1.86 


2.70 + 0.37 


< 0.16 


2359-02 


15 


2.0950 


-0.73 


5.00 + 1.16 


3.94 + 0.26 


0149+33 


15 


2.1400 


-1.64 


3.16 + 0.73 


0.31 + 0.08 


0528-250 


21 


2.1410 


-1.45 


8.91 + 1.03 


1.35 + 0.45 


2359-02 


15 


2.1540 


< -1.03 


2.00 + 0.46 


< 0.80 
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TABLE 10. Summciry of Zn Measurements in DLAs (cont'd) 



QSO 


Ref.t 




[Zn/H] 


iVHi/1020(cm-2) 


iVznii/10i2(cm-2) 


1451+123 


22 


2.2550 


1 


uo 


2.00 ± 0.70 


0.71 ± 0.18 


2348-147 


15,20 


2.2789 






3.63 ± 0.63 


< 0.82 


0216+080 


23 


2.2931 


< 


— n 97 


2 82 ± 1 06 


< 6.46 


PHL957 


15 


2.3091 




"14 


25.10 ± 2.90 


3.12 ± 0.17 


2243-6031 


24 


2.3300 


_i 

X 


OS 

UO 


4.68 ± 0.22 


1.66 ± 0.12 


1232+0815 


25 


2.3376 


n 

u 


S4 


8.00 ± 1.00 


4.90 ± 1.60 


0841 + 129 


15 


2.3745 


1 

- 1 




8 91 ± 1 80 


1 30 ± 15 

_L • V./ vy — 1 — vy • 'Lj 


0112+029 


20 


2.4227 


1 

L 


xo 


9 00 ± 2 00 

ty • vy vy — 1 — • vy vy 


2.86 ± 0.82 


2343+12 


26 


2.4310 


n 

-u 


7/1 


2.24 ± 0.26 


1.73 ± 0.05 


0836+113 


13 


2.4651 


< 


— 1 OQ 


3 80 ± 88 


< 1.32 


1223+178 


13 


2.4658 


_1 




31.62 ± 7.35 


3.55 ± 0.21 


CTQ247A 


27 


2.5505 


_1 




13.49 ± 3.13 


2.75 ± 0.32 


1209+093 


28 


2.5840 


_1 

X 


0^ 


25.12 + 5.84 


9 59 ± 1 17 


CTQ247B 


27 


2.5950 


1 

1 


(14 


12.30 ± 2.86 


4.79 ± 0.22 


2348-011 


13 


2.6145 


< 




19.95 ± 4.63 


< 0.74 


0913+072 


12 


2.6183 


< 


— 1 DQ 


2 30 ± 40 

• vy — 1 — vy • -L vy 


< 0.80 


1759+75 


29 


2.6250 


1 

- 1 


uo 


5.77 ± 0.09 


2.20 ± 0.42 


0812+32 


28 


2.6260 


n 

-u 




22.40 ± 5.20 


10 80 ± 55 

_j_ vy • vy vy — i — vy • v-/ v-/ 


CtQ460 


28 


2.7770 




— 1 97 


10.00 ± 2.32 


< 2.29 


0056+014 


12 


2.7771 


_1 

X 


90 


13 00 ± 2 00 


3 50 ± 2 00 

\j ■ \j vy — 1 — ■ vy vy 


1253-025 


28 


2.7830 


_1 

X 


71 

/ X 


70 79 ± 33 77 


5 86 ± 1 00 

\j ■ v_/ vy — 1 — -1- ■ vy vy 


1132+2243 


28 


2.7830 




— 1 fi4 

— x.u^ 


10.00 ± 1.62 


< 0.97 


1337+113 


28 


2.7950 


< 


X . oo 


8 91 ± 2 07 

vy • ty — 1 — ^ * vy 1 


< 1.59 


2138-444 


8 


2.8520 


1 

- 1 




6.31 ± 1.17 


0.89 ± 0.21 


2342+34 


28 


2 9080 


1 

- 1 


9*^ 


12 59 ± 29 


3 17 ± 80 


1021+3001 


28 


2.9490 




1 1 zL 


5 01 ± 1 16 

\J • vy J- 1 -1- • -1_ \J 


< 1.55 


0426-2202 


28 


2 9830 




— 1 Q7 


31 62 ± 11 14 


< 1.46 


0347-383 


30 


3 0250 


n 


Qfi 


3 63 ± 42 


1 70 ± 20 

_i_ ■ 1 vy — 1 — vy ■ vy 


2334-09 


28 


3.0570 


< 


-0.91 


2.82 ± 0.65 


< 1.47 


0808+52 


28 


3.1130 


< 


-1.15 


4.47 ± 0.72 


< 1.36 


1535+2943 


31 


3.2020 


< 


-0.75 


4.47 + 1.57 


< 3.40 


2344+0342 


28 


3.2190 


< 


-1.71 


22.39 + 3.62 


< 1.85 


2315+0921 


31 


3.2190 


< 


-2.03 


22.39 ± 7.89 


< 0.89 


1506+522 


28 


3.2240 


< 


-1.20 


4.68 + 0.76 


< 1.27 


1432+39 


28 


3.2720 


< 


-1.23 


17.78 + 4.13 


< 4.42 


0957+33 


28 


3.2800 


< 


-0.95 


2.82 + 0.52 


< 1.35 


2239-386 


12 


3.2810 


< 


-1.01 


5.75 + 1.01 


< 2.40 


2155+1358 


28 


3.3160 


< 


-1.13 


3.55 ± 1.25 


< 1.13 
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TABLE 10. Summciry of Zn Measurements in DLAs (cont'd) 



QSO 


Ref.t 


^ahs 


[Zn/H 




ArHi/1020(cm-2) 


iVz„n/10i2(cm-2) 


1715+3809 


31 


3.3410 


< -1.57 


11.20 + 3.47 


< 1.30 


1117-1329 


32 


3.3504 


-1.21 




6.92 + 1.70 


1.81 + 0.12 


0000-263 


33 


3.3901 


-2.03 




25.70 + 4.76 


1.02 + 0.12 


1802+5616 


31 


3.3910 


< -0.53 


2.00 + 0.46 


< 2.54 


0007+2417 


31 


3.4960 


< -1.34 


12.60 + 2.92 


< 2.45 


1802+5616 


31 


3.5540 


< -0.50 


3.16 + 0.74 


< 4.29 


1723+2243 


28 


3.6950 


-0.62 




3.16 + 1.11 


3.22 + 1.07 


1248+31 


28 


3.6960 


< -0.65 


4.27 + 0.69 


< 4.12 


1535+2943 


31 


3.7610 


< 0.08 




2.51 + 0.89 


< 12.90 


0133+0400 


28 


3.7740 


< -0.08 


3.55 + 0.98 


< 12.56 


0747+2739 


28 


3.9000 


< -0.73 


3.16 + 0.74 


< 2.50 



t References: 1. This Paper; 2. Boisse et al. (1998); 3. dc la Varga et al. (2000); 4. 
Meyer & York (1992); 5. Khare et al. (2004); 6. Rao, Turnshck, & Nestor (2004); 7. Pettini 
et al. (2000); 8. Prochaska et al. (2003b); 9. Meyer et al. (1995); 10. Pettini et al. (1999); 
11. Lopez et al. (1999); 12. Pettini et al. (1997); 13. Prochaska et al. (2001b); 14. Ellison 
& Lopez (2001); 15. Prochaska & Wolfe (1999); 16. Petitjean, Srianand, & Ledoux et al. 

(2000) ; 17. Ledoux, Srianand, & Petitjean (2002); 18. Petitjean, Srianand, & Ledoux et al. 

(2002) ; 19. Prochaska & Wolfe (1998); 20. Pettini et al. (1994); 21. Centurion et al. (2003); 
22. Dessauges-Zavadsky et al. (2003); 23. Lu et al. (1996); 24. Lopez et al. (2002); 25. Ge, 
Bechtold, & Kulkarni (2001); 26. Dessauges-Zavadsky et al. (2004); 27. Lopez & EUison 

(2003) ; 28. Prochaska et al. (2003c); 29. Prochaska et al. (2002); 30. Ledoux, Petitjean, & 
Srianand (2003); 31. Prochaska et al. (2003a); 32. Peroux et al. (2002); 33. Molaro et al. 

(2001) . 
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TABLE 11 

Global Zn Metallicity vs. Redshift (Binned) 



z range 


Med. 

z 


T t 


log (Z/Zq) 
Max. Lim. 
Zn only 


1 f ry 1 ry \ 

log {Z/Zq) 
Min. Lim. 
Zn only 


1 ( ry 1 ry \ 

log (Z/Zq) 
Min. Lim. 
Zn±others 


1 ( ry 1 ry \ 

log (Z/Zq) 
Mean of 
Max, Min'f'f 


1 ( ry 1 ry \ 

log {ZjZc^ 
Surv. Anal. 


0.09-1.37 


0.74 


10, 5 


-0.86 


-1.01 


-0.93 


-0.89 


-0.94 








±0.11 


±0.14 


±0.10 


±0.12 


±0.16 


1.42-2.00 


1.88 


14, 


-1.07 


-1.07 


-1.07 


-1.07 


-1.07 








±0.17 


±0.17 


±0.17 


±0.17 


±0.16 


2.02-2.34 


2.14 


10, 5 


-1.09 


-1.19 


-1.14 


-1.12 


-1.16 








±0.13 


±0.10 


±0.10 


±0.12 


±0.23 


2.37-2.78 


2.59 


10, 4 


-1.20 


-1.25 


-1.22 


-1.21 


-1.21 








±0.08 


±0.10 


±0.08 


±0.08 


±0.14 


2.78-3.27 


3.02 


4, 11 


-1.48 


-1.92 


-1.69 


-1.59 


-1.65 








±0.14 


±0.23 


±0.14 


±0.20 


±0.19 


3.28-3.90 


3.44 


3, 11 


-0.86 


-1.80 


-1.43 


-1.14 


-1.30 








±0.25 


±0.51 


±0.22 


±0.46 


±0.15 



f D: Number of Zn Detections, L: Number of Zn Limits, ff: The minimum limit 
value used while calculating this mean is the one with information from other elements used 
in cases of Zn limits. 
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TABLE 12 

Exponential Fits to Global Metallicity-Redshift Relation 



Redshift range 


Binning 


Treatment of Limits 


log(Zo/Z0) 


B 


0.09-3.90 


Binned 


Maximum limits 


-0.74 ±0.15 


-0.18 ±0.06 


0.09-3.90 


Binned 


Minimum limits, Zn only 


-0.75 ±0.18 


-0.22 ±0.08 


0.09-3.90 


Binned 


Minimum limits, Zn -|- others^ 


-0.71 ±0.13 


-0.23 ±0.06 


0.09-3.90 


Binned 


Mean of Max., Min. Limits^ 


-0.72 ±0.16 


-0.20 ±0.07 


0.09-3.90 


Binned 


K-M Survival Analysis 


-0.79 ±0.18 


-0.18 ±0.07 


0.09-3.28 


Binned 


Maximum limits 


-0.66 ±0.15 


-0.23 ±0.07 


0.09-3.28 


Binned 


Minimum limits, Zn only 


-0.77 ±0.18 


-0.22 ±0.08 


0.09-3.28 


Binned 


Minimum limits, Zn -|- others^ 


-0.69 ±0.14 


-0.24 ± 0.06 


0.09-3.90 


Unbinned 


Maximum limits 


-0.86 ±0.15 


-0.11 ±0.08 


0.09-3.90 


Unbinned 


Minimum limits, Zn only 


-0.95 ±0.08 


-0.14 ±0.05 


0.09-3.90 


Unbinned 


Minimum limits, Zn -|- others ^ 


-0.86 ± 0.08 


-0.15 ±0.05 



f : Here, the minimum limit calculation refers to the one with information from other 
elements used in cases of Zn limits. 



